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The quantum version of the Bochkov-Kuzovlev identity is derived on the basis of the 
appropriate definition of work as the difference of the measured internal energies 
of a quantum system at the beginning and at the end of an external action on 
the system given by a prescribed protocol. According to the spirit of the original 
Bochkov-Kuzovlev approach, we adopt the "exclusive" viewpoint, meaning that the 
coupling to the external work-source is not counted as part of the internal energy. 
The corresponding canonical and microcanonical quantum fluctuation theorems arc 
derived as well, and are compared to the respective theorems obtained within the 
"inclusive" approach. The relations between the quantum inclusive-work w, the 
exclusive- work wq and the dissipated- work Wdis, are discussed and clarified. We 
show by an explicit example that wq and Wdis are distinct stochastic quantities 
obeying different statistics. 
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1. Introduction 

One of the main objectives of nonequilibrium thermodynamics is the study of the 
response of physical systems to applied external perturbations. Around the mid 
of the last century major advancements were obtained in this field with the de- 
velopment of linear response th e ory b y severa l auth ors, among which we mention 
Callen fe Welter] ill 95 lh: iGreenl (|l952h : iKubol (Il957l) . This th eory, insp i red b y the 
works of lEinsteir] (|l926l ) on the Brownian movement and of Johnson ( 19281 ) and 
Nvauist ( 19281) on noise in electrical circuits, established that, under certain circum- 
stances, the linear response to small perturbations is determined by the equilibrium 
fluctuations of the system. In particular, the linear response coeffici ents are prop or- 
tional to two-point correlation functions for Hamiltonia n systems (|Kubo L 19571) as 



well as for stochastic, generally non-equilibrium systems ( Hanggi fe Thomas! 19821) . 
In principle, an infinite hierarchy of higher order fluctuation-dissipation relations 
connects the n-th order response coefficients to (n + l)-point correlation functions. 

In contrast, fluctuation theorems are compact relations that provide information 
about the fully nonlinear response. A ccordingly, fluctuation-dissipation r elations of 



all orders can be derived therefrom. iBochkov fc Kuzovlevl ( 19771 11981 ) were the 



first to put forward one such fully nonlinear fluctuation theorem. These authors 
noticed that, for a classical system, their general fluctuation theorem implies the 
following, extremely simple nonequilibrium identity: 



-PW 



> = 1 



(1.1) 



Article submitted to Royal Society 



TgX Paper 



2 



M. Campisi, P. Talkner and P. Hdnggi 



where Wq is the work done on the system by the external perturbation during 
one specific realization thereof, (•) denotes the average over many realizations of 
the same perturbation, and f3 — (fc^T) -1 , with T the initial temperature of the 
system and ks Boltzmann's constant. Due to the properties of convexity of the 
exponential function, an almost immediate consequence of (jl.ll) . is the second law 
of thermodynamics in the form, (Wo) > 0; i.e. when a system is perturbed from an 
initial thermal e quilibrium, on average , it can only absorb energy. 

The works of lBochkov fc Kuzovlevl ( 19771 Il981 ) has recently re-gained a great 



deal of attention, after Jarzvnskil (jl997fl derived, within the framework of classical 
mechanics, a salient result similar to Eq. (jl.lj) 



-pw\ - 



-PAF 



(1.2) 



which, in contrast to Eq. (jl.l[) , allows to extract an equilibrium property of the sys- 
tem, i.e., its free energy (difference) F, from nonequilib r ium f luctu ations of work W. 
Evidently, th e definitions of work adopted bv l Jarzvnskil (|l997 ) and Bochkov fe Kuzovlev 
( 1977 . 198lh (denoted here respectively as W and Wo) do not coincide. The rela- 
tionships between these two work definitions and the corresponding nonequilibrium 
i dentities, Eas. U.l] ! .2 j. were discussed in a very clear and elncidai in;; manner in 



Jarzvnski fc Horowitzl 2007 ). which, for the sake of clarity, we 



([Jarzvnskil 1200 
summarize below. 

Let us express the time dependent Hamiltonian of the driven classical system 
as the sum of the unperturbed system Hamiltonian Hq and the interaction energy 
stemming from the coupling of the external time dependent perturbation X(t) to 
a certain system observable Q: 



H(q, p; t) = H (q, p) - X(t)Q(q, p), 



(1.3) 



We restrict ourselves to the simplest case of a protocol governed by a single "field" 
X(t) coupling to the conjugated generalized coordinate Q. Generalization to several 
fields Xi(t) coupling to different generaliz ed coordinates Q j is straightforward. 

The definition of work W according to I Jarzvnski (1997) stems from an inclusive 
viewpoint, where one counts the term X{t)Q as being part of the system internal 
energy . In c ontrast the definition of work Wo according to iBochkov fc Kuzovlev 
(|1977l . 119811) belongs to an exclusive viewpoint where instead, such term is not 
counted as part of the energy of the system. More explicitly, if qo , po is a certain 
initial condition that evolves to q/, p/ in a time tf—to, according to the Hamiltonian 
evolution generated by H , then the two different definitions of work become: 

W = ff(q / ,p / ;t/)-.ff(qo,Po;io) (1-4) 
W = H (q f ,p f )-H (q 0)Po ) (1-5) 

It is important to stress that IBochkov fc Kuzovlev (Il977l Il98lh only obtained 
Eq. (Il.lj) in the classical case, notwithstanding the fact that they developed a quan- 
tum version of t heir theory, as well. This diff i culty is related to the fact that work 
was identified by IBochkov fc Kuzovlevl (|l977t 1 1 9 8 lh with the quantum expectation 
of a pretended work operator, given by the difference of final and initial Hamiltonian 
in the Heisenberg representation. To be more clear, if the quantum Hamiltonian 
reads: 

H{t) =H - X(t)Q (1.6) 
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where now H, Hq and Q are hermitian operators, the work operator was defined by 
Bochkov fe Kuzovlevl (|l977l Il98ll) as: Q 



W = H^{tj)-H a 



(1.7) 



where the superscript H denotes Heisenberg picture. A simila r approach was em- 

ploye d also within the inclusive viewpoint, with work defined as (jAllahverdvan fc Nieuwenhuizenl 
20051) 

W = H H (t f )-H (1.8) 

As pointed out clearly by some of us before with the work in (jTalkner et all . 120071 ) 
the Jarzynski Equality (jl.2D cannot be obtained on the basis of the work operator 
(|1.8p . Likewise the Bochkov-Kuzovlev identity ([l.ip cannot be obtained on the 
basis of Eq. (|1.7I) . It is by now clear that the impossibility of extending the classical 
results p. 111. 21) on the basis of quantum work operators (|1.7ll.8p . respectively, is 
related to the fact that work characterizes a process, rather than a state of the 
system, and consequently cannot be represented by an observable that would yield 
work as the result of a single projective measurement. In contrast, energy must be 
measured twice in order to determine work, once before the protocol starts and a 
second time immediately after it has ended. The difference of the outcomes of these 
two measur ements then y i elds t he work performed on the system in a particular 



realization (jTalkner et at 
In this 



2007) 



(1977 



aper we will adopt the exclusive viewpoint of iBochkov fe Kuzovlev 
19811 ). but use the proper definition of work as the difference between the 



outcomes of two projective measurements of Hq, to obtain the quantum version of 
Eq. (|l.ip . Indeed we will develop the theory of quantum work fluctuations within the 
exclusive two-point measurements viewpoint in great generality. In Sec. [2] we study 
the characteristic function of work. In Sec. [3] and Sec. IHw e derive the exclusive ver- 
sions of the Tasaki-Crooks quantum fluctuation theorem (jTalkner fc Hanggil . 12007 : 
Talkner. Campisi fc Hanggi , 120091 : ICampisi et aZ.I. l2009h. and of the microcanoni- 
cal quantum fluctuation theorem ( Talkner. Hanggi fc Morillol 20081 ). respectively. 
Sec. [5] closes the paper with some remarks concerning the relationships between the 
inclusive-work, the exclusive-work, and the dissipated-work. 



2. Characteristic function of work 

As mentioned in the introduction, work is properly defined in quantum mechanics 
as the difference of the energies measured at the beginning and end of the protocol, 
i.e., at times to and f/ > to, respectively. According to the exclusive viewpoint that 
we adopt here, this energy is given by the unperturbed Hamiltonian Hq. Let e n and 
|n, A), denote the eigenvalues and eigenvectors of Hq: 

H \n,X) = e n \n,X). (2.1) 

Here n is the quantum number labeling the eigenvalues of Hq and A denotes further 
quantum numbers needed to specify uniquely the state of the system, in case of 

f IBochkov fc Kuzovlevl (1 1977T) defined the "operator of energy absorbed by the system" E = 
f*£ X(t)Q h (t)cIt, where Q H (r) is the operator Q in the Heisenberg representation. It is not 
difficult to prove that E coincides with Wo in Eq. 111.71 1. 
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degenerate energies. A measurement of Hq at time to gives a certain eigenvalue e n 
while a subsequent measurement of Hq at time tf gives another eigenvalue e m , so 
that 

wo = e m - e„ (2.2) 

Evidently wo is a stochastic variable due to the intrinsic randomness entailed in the 
quantum measurement processes and possibly in the statistical mixture nature of 
the initial preparation. In the following we derive the statistical properties of wq, 
in terms of its probability density function (pdf), and the associated characteristic 
function of work. 

Let the system be prepared at time t < to in a certain state described by the 
density matrix p(to). We further assume that the perturbation X(t) is switched 
on at time to- At the same time the first measurement of Hq is performed, with 
outcome e n . This occurs with probability: 

p n = J2{n,X\p(t )\n,X) =TrP„p(i ) (2.3) 

A 

where P n is the projector onto the eigenspace spanned by the eigenvectors belonging 
to the eigenvalue e„: 

P n = ^hA)<n,A| (2.4) 

A 

and Tr denotes trace over the Hilbert space. According to the postulates of quantum 
mechanics, immediately after the measurement the system is found in the state: 

Pn = P n p{to)P n /p n ■ (2.5) 

For times t > to the system evolves according to 

P(t) = U tMPn ul t0 (2.6) 

with Utj, denoting the unitary time evolution operator obeying the Schrodinger 
equation governed by the full time dependent Hamiltonian (Eq. II .6( 1: 

ihd t U t , t0 = H{t)U t , t0 , U t0 , t0 = 1 . (2.7) 

At time tf the second measurement of Ho is performed, and the eigenvalue e m is 
obtained with probability: 

p(m\n) = TrP m p n (t f ). (2.8) 
Therefore the probability density to observe a certain value of work wq is given by: 
Pt f ,to( w o) = ~ t em ~ e n ])p(m\n)p n . (2.9) 

m,n 

We use the superscript throughout this paper to indicate the exclusive viewpoint. 
The same symbols, without the superscript denote the respective quantities within 
the inclusive viewpoint. 

The Fourier transform of the probability density of work gives the characteristic 
function of work 

G° tfM {u) = J dw p° tftto (w )e™ w ° (2.10) 
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which allows quick derivations of fluctuation theorems an d nonequilibrium equali- 



ties. P erforming calculations analogous to those reported bv lTalkner. Hanggi fc Morillo 



(2008) we find the characteristic function of work, in the form of a two point quan- 



tum correlation function: 

G? ^(u) = Tr^ uH oM e - iuHo p{t ) = (e iuH °^e- iuHo ) (2.11) 
where p(to) is defined as: 

n n 

and the superscript H stands for Heisenberg representation, i.e.: 

<(</) = Ut tM H Q U tf , t0 (2.13) 

This exclusive-work characteristic function G° tg should be compared to the in- 
clusive-work characteristic function Gt s ,t that is obtained when looking at the 
difference w of the outcomes E n (to) and E m (tf) of measu rements of th e tota l 



time dependent Hamiltonian H( t). In this case one obtains ([Talkner et all 12007 



Talkncr. Han ggi fc Morillol 1200 



G tf , t0 (u) = Tr^ uHH We- iuHo p(to) = {e luH " ( t f) e -™ H °) (2.14) 

The difference lies in the distinct fact that (tf) appears in the exclusive approach 
in place of the full H H (t f ). 

(a) Reversed protocol 
Consider next the reversed protocol 

X(t) = X(t f +t -t) (2.15) 

which consecutively assumes values as if time was reversed. Let H (t) be the resulting 
Hamiltonian: 

H(t) =H - X(t)Q (2.16) 
The characteristic function of work now reads: 

G t/ ,to(u) = Tre iuSH ^e- iuHo p(t ) = ( e ™^?(*/) e -»^o) ( 2 .17) 

where 

H"(t f ) = Ul fM H»U ti , t0 (2.18) 
and Ut ft t is the time evolution operator generated by H{t): 

ihd t Ut,to = H{t)U tM U toM = 1 (2.19) 
Assuming that the Hamiltonian H{t) is invariant under time reversal i.e.0 

OH^G- 1 = H{t) , (2.20) 

f Here we assume that the Hamiltonian does not depend on any odd parameter, e.g., a magnetic 
field. Treating that case is str aightforward and amou nts to reverse the sign of the odd parameter 
in the r.h.s. of Eq. lEOOl) . see llAndrieux et aLl.l2009h . 
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where O is the antiunitary time reversal operator ([Messiah! . 119621 ) , the time evolu- 
tion operators associated to the forward and backward protocols are related by the 
following important relation, see |Appcndix A| 



Ut ,t f = U} }tt0 = QU tftt0 ®- x . 



(2.21) 



In the following section we will derive the quantum version of Eq. and its 

associated work fluctuation theorem. This will be accomplished by choosing the 
initial density matrix to be a Gibbs canonical state. In Sec. 2] we will, instead, 
assume an initial microcanonical state. 



3. Canonical initial state 

For a system staying at time to m a canonical Gibbs state: 

p(t ) = p(t ) = e-P Ho /Z (3.1) 

where Zq = Tre - ^ , p(to) coincides with p(to) because the latter is diagonal with 
respect to the eigenbasis of Hq (see Eq. I2.12j) . Plugging Eq. (|3.1I) into (|2.11[) , we 
obtain: 

G° tfAo (/3;u) =TTe mH ^ t ^e" mHo e-l iHo /Z (3.2) 

where for completeness we have listed the dependence upon /? among the arguments 
of Gj tQ . The quantum version of Eq. (|1.1[) immediately follows by setting u — i(3: 

( e -^> = G° tftto (j3;ij3) = Tre- 0H "W/Z o = Tr e - 0Ho /Z o = 1 (3.3) 

where in the third equation we have used Eq. (|2.13p , the cyclic property of the trace 
and the unitarity of the time evolution operator: U} ft t Q Utf,t = 1- 

Moreover we find the following important relation between G° to and G° to , 
see |Appendix B| 

G° tfito (^u)=G° tjito (0-,-u + il3). (3.4) 

By means of inverse Fourier transform, the following quantum Bochkov-Kuzovlev 
fluctuation relation between the forward and backward work probability density 
functions is obtained: 

^. (3.5) 



This must be compared to the quantum Tasaki-Crooks relation that is obtained 
within the inclusive viewpoint (jTalkner fc Hanggil I2007I) : 



Pt f ,toW; -w) 



(3.6) 



where, in contrast to Eq. (1531) the term AF = -fi^ 1 [In Tr e' pHW - lnTre"' 3 ^], 
appears. 
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(a) Remarks 

Eqs. (13.31 13.51) constitute original quan tum results that do not appear in the 
works of iBochkov fe Kuzovlev ( 19771 [l98lh . In the classical case they found a fluc- 
tuation theorem similar to Eq. (|3.5I) . reading: 



P[Q(t);X(t)} 
P[eQ(t);sX(t)} 



exp 



/? / X(t)Q(t) 

to 



(3.7) 



where P[Q(t); X(t)] is the probability density functional to observe a certain tra- 
jectory Q(t) given a certain protocol X(t). Here Q(t) is a short hand notation for 
Q(q(qo,Po,T),p(qo,Po,r)), see Eq. (O, where (q(q , Po, r), p(q , p , r)) is the 
evolved initial condition qo, po at some time r G [to, t/ ], for a certain protocol -X\t). 
The symbol e denotes the parity of the observable Q under time reversal (assumed 
to be equal to 1 in this paper). The symbol ~ denotes quantities referring to the 
reversed protocol. The classical probability of work p$ * t„ (Wo) is obtained from the 
Q-trajectory probability density functional P[Q(t); X(t)] via the formula: 



VQ{t)P[Q{t)-X(t)}8 



W 



X{t)Q(t) 



(3.8) 



where the integration is a functional integration over all possible trajectories such 
that //' X(t)Q{t) = Wo- With this formula one finds from Eq. (|3.7[) the exclusive 
yersion of the classical Crooks fluctuation theorem for the work probability densities 
( Jarzvnski fc Horowitz! . 12007 ) 



~cl,0 (q 

Pt f ,t (fc 



-W )e 



PW 



(3.9) 



Notably, a quantum version of Eq. Q3.7P does not exists because:... "in the quan- 
tum case it is impossible to i ntroduce unambiguously a [...] probability functional" 



(jBochkov fc KuzovlevL 119811) . It is only by giving up the idea of true quantum 



trajectories and embracing instead the two-point measurement approach that the 
quantum exclusive fluctuation theorem Eq. (13.5[) can be obtained, and has been 
obtained here, for the first time. 



4. Microcanonical initial state 

We consider next an initial microcanonical initial state of energy E, that can for- 
mally be expressed as: 

p(t )=p(t ) = S(Ho-E)/n (E), (4.1) 

wherein Q,q(E) — Tr5(Ho — E). Actually one has to replace the singular Dirac 
function 5{x) by a smooth function sharply peaked around x = 0, but with infinite 
support. A normalized gaussian with arbitrarily small width serves this purpose 
well. 

With this choice of initial condition, the characteristic function of work reads: 

G° tfito (E;u) = Tre^^e-^^tfo - E)/Q (E) 

= Tie tu[H o^^-^S(Ho - E)/n Q (E) (4.2) 
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where for completeness we listed the dependence upon E among the arguments of 
G® to . By applying the inverse Fourier transform we obtain: 

p° t/!to (E;w ) = Tr6(H?(t f ) -E- w )S(H a - E)/Q (E) (4.3) 
Likewise, for the reversed protocol, 

$, ito (E;wo) = TrS(H?(t f ) -E- w )8(H - E)/Q (E) (4.4) 

is found. 

We then find the following relation between the forward and backward work 
probability densities, sec jAppcndix C| 

n (E)p° tftto (E; w ) = n (E + wo)$ />to (E + w ; -w Q ) (4.5) 

Then, the quantum microcanonical fluctuation theorem reads, within the exclu- 
sive viewpoint: 

Pt f ,t (E;w ) n a (E + w ) 



jP tf<to (E + w ;- Wo ) n (E) 



(4.6) 



This must be compared to the quant um microcanonical fluctuat i on th eorem, ob- 
tained within the inclusive viewpoint ( Talkner. Hanggi fe Moriilol 2008 ) 



p tf<to (E;w) _Q f (E + w) 



Pt f ,t (E + w; -w) Qo(E) 

The difference lies in the fact that within the exclusive viewpoint the densities of 
states at the final energy E + wq , is determined by the unperturbed Hamiltonian, 
i.e., Q.q(E+wo) = Ti S(Ho — (E+wq)), whereas it results from the total Hamiltonian 
in the inclusive approach: Qf(E + w) = Tr S(H(tf) — E — w). 



Eq. (|4.7I) was first obtained within the classical framework by (|Cleuren et al. 



20061 ) . It is not difficult to see that Eq. (|4.6[) holds classically as well. 



(a) Remarks 

Just as Eq. (13.511, t his Eq. (14.61) is a new result that was not reported before by 



Bochkov fc Kuzovlevl (|1977l II981I ). It is very interesting to notice, however, that 
those authors already put forward a classical fluctu ation theorem for the micro- 
canonical ensemble, which can be recast in the form (jBochkov fc Kuzovlev . 1981 ): 



P[I(t);X(t);E] _ n (E + W q ) 



P[-eT(t);sX(t);E + W } n o(E) 

where P[J(r); X(r); E] is the probability density functional to observe a certain 
trajectory J(r) given a certain protocol and an initial microcanonical ensemble of 
energy E. Here 

J(t) = Q(q(qo, po, r), p(q , Po, r)) (4.9) 

denotes the current. By functional integration the classical microcanonical theorem 
for the pdf of work 

pf; tQ (E,Wo) _n (E + w ) 



pi;° to (E + w ,w ) ME) 
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is obtained from (|4.8I) in the same way as (|3.5I) follows from (13. 7[) . However the quan- 
tum version of (|4.8p does not exists and the derivation of the quantum microcanon- 
ical fluctuation theorem (|4.6p is indeed only possible if the two-point measurement 
approach is adopted. 

The fluctuation relations of Eqs. (|4.6l 14. 7[) can be further exp ressed i n term s of 
entropy, according to the rules of statistical mechanics. Following iGibbsl (|1902l) two 
different prescriptions are found in textbooks to obtain the entropy associated to 
the microcanonical ensemble: 



s(E) = k B In Q(E) = Tr S(H - E) 
S{E) = k B ln$(E) =Tr9(H-E) 



(4-11) 
(4-12) 



The two definitions coincide for large systems with short range interactions among 
their constituents, but may substantially differ if the size of the system under 
study is small. It is by now clear that, of the two, only the second - cus tomarily 
called "Hertz entropy" - is t h e fundamentally correct one (iHertzl (ll910df^:ISchluterl 



Pearson et a/1 (ll985n : ICamp"isil (l2005ll2008il2010l ): lDunkel fc Hilberd (12006^ 



] Using the microcanonical expression for the temperature ksT(E) = (dS(E) / dE) -1 
<&(E) /Q,(E), we can re-express the quantum microcanonical Bochkov-Kuzovlev fluc- 
tuation relation in terms of entropy and temperature as: 



P° tf , t0 (E;w ) _ 

$ ftto (E + w ;-w ) T (E + woY 



ME) 



D [S (E+iu )-S (E)}/k E 



(4.13) 



where the subscript in T and S denotes that these quantities are calculated for 
the unperturbed Hamiltonian Hq. Likewise, adopting the inclusive viewpoint one 
obtains: 



Pt f ,t (E;w) 



ME) 



Pt f ,t (E + w;-w) T f (E + w) 



„[S f (E+w)-S (E)]/k E 



(4.14) 



where the subscript / in T and S denotes that these quantities are calculated for 
the total final Hamiltonian i?(i/).[j] 



5. Discussion 

We derived the quantum Bochkov-Kuzovlev identity as well as the quantum canon 
ical and microcanonical work fluctuation theorems within th e exclusive approach, 



and have el ucidated their relations to the original works of iBochkov fe Kuzovlev 



(|19771 . 119811 ). The extension of the corresponding classical theorems to the quan- 
tum regime is only possible thanks to the proper definition of work as a two-time 
quantum observable. We close with two comments: 

1. For a cyclic process, X(tf) — X(to), inclusive and exclusive work fluctuation 
theorems coincide. However in no way is it true that the exclusive approach of 

f It is interesting to notice that Eins tein was well a ware of the works o f lHertJ il910d rg) which 
he praised as excellent ( "vortrefflich" ) jEiasteirJ. Il91lf) . 

X If instead of the microc anonical ense mble 1 14.111 , the modified microcanonical ensemble g(to ) = 
8(E — Ho)/[TrO(E — Ho)], (Rucllc, 1969) would be used as the initial equilibrium state, then the 



fluctuation theorem assumes the s ame form as in Eq. 114. 1411 . but without the ratio of temperatures 
jTalkner. Hanggi fc Morilld. 1200^1 . 

If Similar remarks were made also within the classical framework by Jarzvnski ( 2007). 
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Bochkov & Kuzo vlev, adopted here, is restricted to cyclic processes, as some authors 
have suggested (|Allahverdyan fc Nieuwenhuizen . 2005 ; Cohen fc Mauzerall , 2005 ; 
Andrieux fe Gaspardl . 120081) . As stressed in the introduction, the difference of the 



two approaches originates from the different definitions of work, and is not related 
to whether the process under study is cyclic or is not cyclic. 

2. Within the inclusive app r oach i t is natural to define the diss i pated work as 
Wdis — w — AF ( Kawai et all l2007t IVaikuntanathan fc Jarzvnski , 2009h . Then, 
the Jarzynski equality (jl.2l) can be rewritten as (e~^ Wdis ) = 1. This might make 
one believe that the exclusive work wq coincides with the dissipated work Wdis- 
This, though, would be generally wrong. The dissipated work Wdis is a stochastic 
quantity whose statistics, given by p^ s tQ (w dts ) = p t} ,t ( w <Hs + AF), in general 
does not coincide with the statistics of exclusive work wq, given by Pt fl t { w o)- Sec 
|Appcndix D| for a counterexample. Only for a cyclic process, for which AF = 0, 
does the dissipated-work Wdis coincide with the inclusive-work w, which in turn 
coincides with the exclusive- work wq. 
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Appendix A. Derivation of Eq. ( I2.21j) 

The time evolu tion operator Ut,t can be expressed as a time ordered product 
(|Schleichl . [200l : 



U tftn = lim e -mt»)T e -iH(t N - 2 )T __ e -iH(t!) 7 



(Al) 



where t — (tf — to)/N, and t v = to + vt : for v = 0, . . . , N (note that tjy — tf). Due 
to Eqs. ([2~T5l[2~T6)l . it is H(t) = H(t f + t - t), then: 

lim e-iH^i^e-iH^r _ _ _ e -iH(t N )r 



Ui 



tf,to 



(A 2) 



Therefore: 

eutftoQ' 1 ^ lim ee-^^e^ee-^^e -1 . . . ee-^^e -1 (A3) 



where we inserted Q 1 = 1, N — 1 times. Due to the property (|2.20[) and the 
antilinearity of O it is 

0e -iff(t)r e -i = e iH(t)T _ ( A4 ) 

Using this equation, we find: 



eu tttn e~ 1 = iim e iH(ti)T e iH(t 2 )r ^ e w(t N )T 



lim 



N- 



e -iH{t N )r _ _ _ e -iH(t 2 )r e -iH{t!)i 



U t f ,t ~ U t ,t, ■ 



(A 5) 
(A 6) 

(A 7) 
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In a similar way we also obtain: 

u tf ,to = e^e- 1 . (A 8) 

Appendix B. Derivation of Eq. (13. 4p 

The exclusive- work characteristic function reads (|2.11l) 

G° /it0 G8;«) = TrU tf>to e iuH °Ut ftto e- iuB <>e-l >H °/Z . (B 1) 

Using Eqs. (|A 7IA 81) . we obtain: 

G° tftto (p;u) = Treul^Q^e^QUt^Q^e-^e-^Qe^/Zo 

where we have inserted 60 _1 = 1 at the right end. By multiplying Eq. (|A 4[) by 
G _1 to the left and by O to the right, we have (replacing r with u) 

Q-l e iH{t)uQ = e -iH(t)u ) ( B2 ) 

therefore: 

G° tf>t0 (p;u) = TrGUl jto e- iuH °U tf , to e iuH °e-P H °Q- l /Z . (B3) 
The antilinearity of implies, for any trace class operator A: 

TrOAO- 1 = TrA* . (B4) 

Therefore: 

G° tf>t0 (p; u) = Tre-P H °e- mH °ul to e mH °U tf , t0 /Z (B 5) 

Using the cyclic property of the trace finally leads to 

G° tftto (0;u) = IV^^e'C-^^^t jtoe -i(-t^«Ho e -^Ho /Zo (B6) 
= G° /Aj (/3;- U + */?). (B7) 

Appendix C. Derivation of Eq. ( 14.51) 

The microcanonical exclusive- work probability density function reads, Eq. (|4.3p : 

p° ttM {E; w ) = Tr6(H?(t f ) -E- w )6(H - E)/Ci (E) (C 1) 

= TrU} fit0 6(H a -E- w )U tftto S(H - E)/Q (E) (C2) 

Employing Eqs. (I A 7IA 8|) . then leads to: 

Q (E)p° tf to (E;w ) = Tr6LT t/!t0 9-^(^0 -E- Wo)QU^ ^^(Ho - E)^- 1 

(C3) 

where we have inserted 00 _1 = 1 at the end. Being the Dirac delta a real function 
we have 

e- 1 S(H -E)e = S(H -E) (C4) 
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because Hq is assumed to be invariant under time reversal. Then: 

n (E)p° tfjto (E;w ) = TrGU tj , to 6(H -E- w n )u} fAo 5{H - E)^ 1 . (C5) 

Using Eq. (|B 41) . we obtain: 

n (E)p° tf>to (E;w ) = Tr 5{H - E)U tfMl S(H Q — E — w )U} fM . (C6) 

Thanks to the cyclic property of the trace one finally arrives at: 

n (E)p° tf>to (E;w ) = TrU} f>to 6(H - E)U tf , t J(H - E - w ) (C7) 
= n (E + wo)f tjito (E + wo;w ). (C8) 

Appendix D. Comparison between dissipated- work and 

exclusive-work pdf's 

In this appendix we provide an example that shows that the dissipated- work Wdis 
and the inclusive work wq are distinct stochastic quantities with different statisti- 
cal properties. To this end we show that their probability density functions may 
have different first and second moments. We consider a driven quantum harmonic 
oscillator of unit mass and unit angular frequency: 

H(t) =p 2 /2 + q 2 /2- X(t)q (Dl) 

For simplicity we assume to = 0, X(to) = 0, and we chose units in such a way that 
h = 1. Let \n, t) denote the instantaneous eigenvectors of H (t) corresponding to the 
instantaneous eigenvalues E n {t) = (n+ 1/2) — X 2 (t)/2. 

(a) The probability density of dissipated-work 

The probability density function (pdf) of inclusive work, corresponding to an 
initial canonical state, is 



PL 



oH = S ( w - m + n + X 2 (t)/2)\a mn \ 2 e-^ n+1 ^ /Z(Q) (D 2) 



where Z(Q) = ^2 n e ^("+ 1 / 2 ) is the initial partition function, and |a m „| 2 are the 
probabilities to make a transition between two eigenstates of the total Hamiltonian 

\a mn \ 2 = |<m,i|[/ t , |n,0)| 2 (D3) 

where we h ave set to — and tf — t. According to Talkner. Burada fc HanggH 
i;he mean value and the variance of the inclusive work pdf (|D 2|) are 

given by 



(w) = J dxx Pt ,o(x) = L(t) - X 2 (t)/2 (D4) 
(Aw 2 ) = [ dx[x - L(t)] 2 p t ,o(x) = 2UL(t) (D5) 
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where U = J2 n ( n + l/2) e /3( - Tl+1 ^ 2 V^o is the initial average energy, and[j] 

L(t) = C(i) 2 /2 + [S(t) - X(t)} 2 /2 (D 6) 

where 

S(t) = [ dsX(s) sin(* - s) , C(i) = / cZsX(s) cos(i - s) . (D 7) 

Jo Jo 

The partition function of work at the final time t is Z(t) — Z(0)e l3X (*)/ 2 , there- 
fore the free energy difference AF = — \nZ(t)/Z is given by AF = —X 2 (t)/2 
( Talkner, Burada fc Hanggi , 20081) . Hence the dissipated work is 

w dis = w + X 2 (t)/2. (D8) 
Accordingly the dissipated work pdf is 

p%o(w dis ) = p t ,o(w) = p t ,o(wdis - X 2 {t)/2) 

= J2 5 ( w ^ - m + n)\a mn \ 2 e-^ n+1 ^/Z(0) . (D 9) 

mn 

It immediately follows that 

(w dls ) = J dxx p^ (x) = L(t) (D 10) 

(Aw 2 dls ) = J dx[x - L(t)} 2 p^(x) = 2UL(t) . (D 11) 

Note that, as it should be, (wdis) > 0. 

(b) The probability density of exclusive-work 
The exclusive-work pdf is given by 

P°o(«*0 = E - m + n)\a° mn \ 2 e-^ n+1 ^/Z(0) (D 12) 

mn 

where |a^ n | 2 denotes the probability to make a transition between two states of 
the unperturbed Hamiltonian: 



,0 |2 



|(m,0|C/ t! o|n,0)| 2 . (D 13) 



It is known ( Husimi . 19531 : Campisi . 2008^ ) that the transition probabilities |a mn | 2 
depend on the time t at which the second measurement is performed, via the func- 
tion L(t), that is the |a mn | 2 are of the form |a mn | 2 = f nm [L(t)], for certain functions 
f nm that need not be sp ecified here. Usin g Wigner functions to calculate the tran- 
sition probabilities as in ( Campis i. 2008 i, Appendix), we notice that the transition 

f In dTalkner, Burada fc Hanggil, |200S| , l200Sfl L is given as L(t) = | /„' dsf(s)e ia \ 2 , where / = 
—X/y/2. It is a matter of elementary calculus to check that this expression coincides with Eq. 

(nui- 



Article submitted to Royal Society 



14 



M. Campisi, P. Talkner and P. Hdnggi 



probabilities |a^ ln | 2 are obtained from the same expression as of |a m „| 2 , with the 
only difference that Lit) is replaced by 

L (t) = C 2 (t)/2 + S 2 (t)/2 (D14) 

that is |a^„| 2 = f nm [Lo{t)]. Therefore the exclusive- work pdf (|D 12[) is obtained 
from the dissipated- work pdf (ID 9p simply by replacing L(t) with L (t). It follows 
immediately that 

(w ) = J dxxp tt0 (x)=L Q (t) (D15) 

(Aw 2 ) = J dx[x - L {f)] 2 p%{x) = 2UL Q (t) (D 16) 

Note that, as expected, (wo) > 0. 
For the specific protocol 

= 2 sin(i) (D17) 

we find 



L{t) - L (t) = t sin(2£) (D 18) 

which is apparently different from zero except for integer multiples of it/ 2. Thus for 
any duration t of the protocol (ID 17|) that is not an integer multiple of n/2, Lq ^ L. 
Accordingly the first and second moments of pf™ andp° differ, meaning that Wdi 8 
and wq are distinct stochastic variables with different statistical properties. 

It should be stressed that analogous results are found also for a classical driven 
harmonic oscillator. The statistics of dissipated-work and of exclusive-work gener- 
ally differ, this fact holds true both quantum-mechanically and classically. 
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